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Abstract 

Recently there has been renewed interest in second sound in superfluid Bose and Fermi gases. 
By using two-fluid hydrodynamic theory, we review the density response Xnn(q, i^) of these systems 
as a tool to identify second sound in experiments based on density probes. Our work generalizes 
the well-known studies of the dynamic structure factor S'(q, cj) in superfluid ^He in the critical 
region. We show that, in the unitary limit of uniform superfluid Fermi gases, the relative weight 
of second vs. flrst sound in the compressibility sum rule is given by the Landau-Placzek ratio 
eLP = (cp — Cv)/cv for all temperatures below Tc- In contrast to superfluid ^He, clp is much larger 
in strongly interacting Fermi gases, being already of order unity for T ^ 0.8Tc, thereby providing 
promising opportunities to excite second sound with density probes. The relative weights of first 
and second sound are quite different in ^(q, oj) (measured in pulse propagation studies) as compared 
to Imx„„(q,a;) (measured in two-photon Bragg scattering). We show that first and second sound 
in S'(q, w) in a strongly interacting Bose-condensed gas are similar to those in a Fermi gas at 
unitarity. However, in a weakly interacting Bose gas, first and second sound are mainly uncoupled 
oscillations of the thermal cloud and condensate, respectively, and second sound has most of the 
spectral weight in 5'(q, oj). We also discuss the behaviour of the superfluid and normal fluid velocity 
fields involved in first and second sound. 

PACS numbers: 03.75.Kk, 03.75.Ss, 67.25.D- 
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I. INTRODUCTION 

The most dramatic effects related to superffuidity in liquid ^He arise [l| when the dynam- 
ics of the two components are described by the two-fluid hydrodynamics first discussed by 
Landau j2|. These equations only describe the dynamics when the non-equilibrium states 
are in local hydrodynamic equilibrium js], which requires short collision times between the 
excitations forming the normal fluid. (This requirement is usually summarized as ojt ^ 1, 
where u is the frequency of a collective mode and r is the appropriate relaxation rate.) The 
study of ultracold gases when they are in local equilibrium has been difficult because the 
density and the s-wave scattering length are typically not large enough. However, recent ex- 
perimental work on trapped Bose-condensed gases has reported some success, with evidence 
for a second sound mode in highly elongated (cigar-shaped) traps W- 

Another approach to achieving conditions where the Landau two-fluid description is cor- 

n 

rect has been to consider a Fermi superfluid gas close to unitarity [5!], where the s-wave 
scattering length between Fermi atoms in two different hyperfine states is infinite. Develop- 
ing earlier theoretical studies jo, 7|, Taylor and co-workers js] have recently given detailed 
predictions for the first and second sound breathing oscillations in a trapped Fermi gas at 
unitarity. The coupled differential equations of the Landau two-fluid description were solved 
variationally, with results which agreed with analytic predictions at T — and T — )■ T^. So 
far, only an isotropic trap has been studied. 

In jsjl, it was shown that (as in superfluid "^He) the frequencies of first and second sound 
in unitary Fermi gases were quite well approximated by assuming that the solutions of the 
two-fluid equations corresponded to pure uncoupled density and temperature waves. This 
■s ■„ sharp contrast to the situation ■„ dilute, weakly intetaetiug Bose gases, where first 
and second sound involve both density and temperature oscillations [9y. Naively, the first 
result above would lead one to believe that second sound in a unitary Fermi gas would have 
very small weight in the density response function. However, Arahata and Nikuni [l^ have 
shown that in a uniform Fermi gas at unitarity, a density disturbance can have first and 
second sound pulses of comparable magnitude at temperatures of order 0.8Tc. 

The purpose of this paper is to provide a systematic study of the density response func- 



tion x„„(q, cj) of a uniform superfluid atomic gas in the hydrodynamic regime ll|, Il2| . 



as described by the non-dissipative Landau two-fluid equations (see sections [TTl and IHip . 
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In this hydrodynamic region, the related dynamic structure factor is given by 5'(q, u) oc 
Imx„„(q, a;)/ci;. In studies of superfluid "^He, S{q,u) is directly measured in Brillouin light 
scattering experiments Q, Q. In dilute gases, density pulse propagation experiments ef- 



fectively measure S{q,u) as noted in [10|. In addition, experiments using two-photon Bragg 
scattering have a cross-section proportional to Imx„„(q, cj). While all these experiments 
probe the same density response function, Imx„„(q, w), the relative weight of the low fre- 
quency resonances are strongest in ^(q, w) because of the extra factor of 1/u noted above. 
This fact has significant implications for measuring second sound since the frequency of 
second sound becomes much smaller than first sound as is approached. 

We relate our analysis to classic discussions of the dynamic structure factor ^(q, u) in 



13l-ll7l|. However, many new 



the two-fluid region of superfluid ^He in the critical region 
features arise when dealing with quantum gases. 

New aspects of our work include the use of the compressibility sum rule (sections UTT] 
and IIVI) in understanding the relative weights of first and second sound in the density 
response function. We also discuss the role played by the Landau-Placzek ratio eLp = (cp — 
Cv)/cy in determining these weights, with emphasis on new aspects that arise in superfluid 
gases. Here, Cp and are the specific heats per unit mass at constant pressure and volume, 
respectively. We discuss the differences between S{q, u) of a unitary Fermi gas and that 
of a Bose-condensed gas (sections IIVI and IVT) . and present results for Bragg scattering with 
localized beams applied to the centre of a trapped unitary Fermi gas using the local density 
approximation (section IVll) . In [HI we provide a detailed analysis of the superfluid and 
normal fluid velocity fields associated with first and second sound. 



II. DENSITY RESPONSE FUNCTION 



In this section, we review some basic properties of the density response function Xrm(q, i^), 
the related dynamic structure factor ^(q, u), and the experiments that measure these quan- 
tities. The density response function involves the correlation between the density at two 
different points at different times: 

Xnn{r,t) = -zeit){[p{r,t),p{0,0)]). (1) 
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Here, 6{t) is the step function and we have set the volume equal to unity. The Fourier 
transform ImXnn(q, i^) of the imaginary part of the density response function is related to 
the dynamic structure factor S'(q, u) by 

ImXnn(q,w + ^0+) = -nn[S{ci,uj) - S{-q,-uj)] 

= -n7r5(q,u;)(l -e-^'^), (2) 

where n is the density, /3 = 1/T is the inverse temperature, and in the last step we have 
made use of the "detailed balance" relation, S{q, —u) = e~^'^S'(q, w) and also the inversion 
symmetry of the system which gives S'(— q, w) = S{q,u). Here and throughout this paper 
we set h = ks = 1- Equation ([2]) can also be re- written as 

5(q,a;) = - — [iV°(a;) + l]ImXnn(q, w + ^0+), (3) 
nn 

where N^{uj) = {e^'^ — 1)^^ is the Bose distribution function. The latter arises in both 
Bose and Fermi fluids because density fluctuations always obey Bose statistics. We see from 
(j2]) that Imx„„(q, w) is the antisymmetric (frequency) part of the dynamic structure factor 

The dynamic structure factor ^(q, w) is measured by using inelastic scattering (e.g., neu- 
tron and Brillouin light scattering) of probe particles, which transfer momentum q and 
energy u to the system, ^(q, w) is discussed and calculated in all standard texts on many 



body physics (convenient reviews are given in [18|, ll9|]). In contrast, the density response 
function x„„(q, w) describes the density fluctuation induced by an external perturbing po- 
tential. As an example, in atomic gases, the imaginary part of the density response function 
is probed in two-photon Bragg scattering and also describes the density pulse produced by 
a blue-detuned laser. Brillouin and neutron scattering have been used extensively to study 
the collective modes in ^He. These techniques cannot be applied to dilute atomic gases, 
however, since the gases are far too dilute to generate an appreciable signal. Instead, spec- 
troscopic probes of dilute gases have successfully used Bragg scattering since this technique 
makes use of the stimulated light scattering processes induced by two counter-propagating 
laser beams, resulting in a strong enhancement of the signal. 



In Bragg scattering experiments (see, e.g.. 



20j), by measuring the total momentum trans- 



ferred to the sample, one measures the imaginary part of the density response function 
X„„(q, w), where q and uj are the momentum and energy transferred by the stimulated ab- 
sorption and emission of the photons ^]. If the Bragg pulse is short compared to l-n juj^^ 



Method 


System 


Quantity measured 


Response function 


Bragg Spectroscopy 


Gases 


Momentum transferred 


Imx„„(q,tj) 


Neutron Spectroscopy 


^He 


Number of scattered neutrons 


5(q,^) 


Brillouin Spectroscopy (cj <^ T) 


^He 


Number of scattered photons 


S{ci,uj) ~ ^Imx„n(q,w) 


Density pulse 


Gases 


Pulse amplitude 


^Imxnn(q,t^) 



TABLE I: Summary of experimental probes. From left to right: the experimental method, the 
system(s) to which the method can be applied, the quantity being measured, and the density 
response function involved. In the two-fluid region of strongly interacting Fermi gases and ^He, 
second sound is weakly coupled to density fluctuations but has a small frequency uj at high temper- 
atures (below Tc). This means that, as a result of the extra factor of 1/uj multiplying the density 
response function, Brillouin spectroscopy and density pulses are more sensitive to second sound 
than Bragg and neutron scattering. 



where is the frequency of the harmonic trap along the axis of light propagation (q = zq), 
the momentum transferred is 



V 



(4) 



where V is the strength of the potential induced by the lasers and r is the pulse duration. 

As noted recently by Arahata and Nikuni [lO|, the density response function can also be 
measured by exciting a density pulse in a uniform gas. Within linear response theory, the 
density fluctuation 6n{r, t) induced by the application of an external perturbing potential 
6V{r,t) is 

Snir,t) = ^Xnn{^,co)SVici,co)e^^-'^~^-\ (5) 

where SV{ci, u) is the Fourier transform of the perturbing potential. Similar to the sound 
propagation experiments in 



2ll . |22| , consider a localized potential applied for a short du- 
ration r < t < 0, and turned off at t = 0. Assuming for simplicity that the perturbing 



potential only varies spatially along the z-axis (and q = zg), for t > 0, one obtains 



6n{z,t) = — dq f/a;5y(g)_^!^!^^:i^e*9^-*-*. (6) 
In-^ J J CO 

Equation ([6]) reveals an important feature: in contrast to Bragg scattering, where the re- 
sponse is proportional to Imx„„(q, w), the response involved in the excitation of density 
pulses in proportional to Imx„„(q, a;)/^. As we discuss in section UTTl for uniform super- 
fluids, this factor of l/u leads to an enhancement in the signal of second sound in density 
pulse experiments as compared to its signal in Bragg spectroscopy. In table HI we review 
the quantities measured by the experimental probes of interest in both dilute gases and 
superfiuid ^He. 

Before discussing the form of the density response function in the two-fluid hydrodynamic 
region in section IIIIl we first review below some general properties of the dynamic structure 
factor that will be of use later on in understanding the contributions from first and second 
sound in experimental probes of superfiuid atomic gases. 



Quite generally, S{q,u) satisfies various 
particular interest in the two-fluid domain 



requency moment sum rules. Two are of special 



23| . The /-sum rule (valid for all q) is 



g2 

^a;5(q,a;) = — . (7) 

The compressibility sum rule arises from the exact Kramers-Kronig identity (also valid for 
all q) 

Xnn(q,w = 0)=/ , (8) 

which involves an inverse frequency moment. (Note that x„„(q, w = 0) = ReXnn(q, = 0) 
is purely real.) Using ([2]), one can show that ([HD is equivalent to 



X'. 



,„(q, u = Q) = -2n r duj'^^^^. (9) 

J-OQ ^ 



The usefulness of (IHl) and ([9]) is due to the fact that in the long wavelength limit (q — )■ 0), the 
density response function Xnn(q = 0,a; = 0) describes static thermodynamic fluctuations 



which can be related to the equilibrium isothermal compressibility 23|, |24 | 



limxnn(q,w = 0) = -n — 



(10) 



T 

Thus we arrive at the so-called compressibility sum rule for the dynamic structure factor 

,S(q,u;') 1 dp 



lim / duj 



q^o J_^ oj' 2m dp 



(11) 



where we have introduced the isothermal sound velocity vt- 

The sum rules in ([7]) and f|TT]) are general and may be viewed as two constraints which 
S'(q, must always satisfy. The /—sum rule is used frequently in discussions of inelastic 
neutron scattering as a check on the experimental data. As we discuss in section IIIH the 
compressibility sum rule has special significance when we discuss the two-fluid hydrodynamic 



region, as first emphasized by Nozieres and Pines 23 1. 



III. DYNAMIC STRUCTURE FACTOR IN THE TWO-FLUID REGIME 

The Landau two-fluid equations for an isotropic (such that Xnn(q, i^) = Xnn{(l,<^)) super- 
fluid (Bose or Fermi) describe the dynamics when collisions are strong enough to produce 
local equilibrium. Using the non-dissipative two-fluid equations, one finds the density re- 



sponse function (see Refs. 



11 



12| and page 138 of [3]), 

m [iij^ — u1q'^){u'^ — u\q^) 



Xnn{q,Uj)-—— ^ ,2 _„2^2V (^2) 



where we have defined a new velocity 

^2^T^^. (13) 

Cv PnO 

Here, Sq = So/Nm is the equilibrium entropy 5*0 per unit mass and = T{ds/dT)p is 
the specific heat per unit mass at constant volume. The equilibrium superfluid and normal 
fluid densities are denoted by p^o and pno, respectively. In f[T^ . ui and U2 are the well- 
known exact first and second sound velocities given by Landau's hydrodynamics. For later 
purposes, we note that these velocities satisfy the exact relations 

2,2 2,222222'^ f-l A\ 

where vt is defined in ( ITTj) and the adiabatic sound velocity is defined as 

We note that thermodynamic relations show that, quite generally, the ratio of the 
specific heats can be related to the isothermal and adiabatic sound velocities, 

7-? = |- (16) 



The difference between Vg and vt is a consequence of the finite thermal expansion. 

Equation f lT2|) gives the density response for Bose and Fermi fluids described by the 
non-dissipative Landau two-fluid equations. The two-fluid density response function in the 
presence of dissipation arising from transport coefficients has been given by Hohenberg 
15 1 . (See also the discussion given by Vinen [l^.) For our purposes 



and Martin 
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in this paper, the primary effect of dissipation is to broaden the delta function peaks in 
the imaginary part of the non-dissipative density response function. It does not otherwise 
change the structure of the two-fluid response function and in particular, does not affect in 
a significant way the weights of first and second sound in this function. 

From f ll2p . it is straightforward to obtain the imaginary part of the density response 
function, which is the experimentally relevant quantity in dilute gases: 

lYaXnn{q,uj) = -n'K—<-^[5{uj-Uiq) + 5{u + Uiq)] 

+ ^[6{u-U2q) + 6{co + U2q)]y (17) 
Here, we have defined the "traditional" amplitudes for first and second sound, namely 

= — 2 2 ' ^2 = o- \^^) 

Ul — U2 Ml — U2 

Using ( !T7|) in (|3]), the two-fluid expression for the dynamic structure factor is found to be 
S{q, uj) = ^[N\uj) + 1] 1^ - u^q) + 5{uj + u^q)] 

+ ^[5{uj~U2q) + 5{u + U2q)]^. (19) 
The two-fluid dynamic structure factor in ffT9l) turns out to satisfy both the /-sum and 



compressibility sum rules, given by ([7]) and (fTTj) [23|. This makes use of the identity N^{uj) + 
1 -|- N^{—uj) -|- 1 = 1, satisfied by the Bose distribution. Substituting f|T9l) into the left-hand 
side of (j7]), one sees that the amplitudes Zi and Z2 describe the weights of first and second 
sound, respectively, in the /-sum rule. That the sum of these two contributions saturate the 
/-sum rule is a consequence of the fact that Zi + Z2 = 1- Similarly, substituting (IT^ into 
the left-hand side of (fTTj) . one sees the relative contributions of first and second sound to 
the isothermal compressibility are given by Zi/u\ and Z^jw^^ respectively. We will denote 
these weights by Wo 

W^ = %, W2 = ^. (20) 



To see that the compressibihty sum rule is saturated by these contributions from first and 
second sound, we note that Zi/u\ + Z2/U2 = l/fy, as can be seen from f|T^ . 

Since we are working in the low frequency two-fluid hydrodynamic region, we can use the 
fact that uj <^T. In this limit, the detailed balance factor in f lT^ can be approximated by 

ArO(u;) + 1 ~ - + i + . . . (21) 

oj z 

The first quantum correction, given by the 1/2 term in fl?I]) in this expansion, is needed to 
satisfy the sum rules previously discussed. This approximation (w ^ T) is often referred to 
as the "classical limit" in textbooks discussing the dynamic structure factor, although in fact 
it describes both superfiuid (T < Tc) and normal fluid (T > Tc) coUisional hydrodynamics. 
Keeping only the leading order first term in (pTj) . ( IT9|) simplifies to 

T ^ Z 

S{q,uj) ^ —<^[5{uj-uiq) + 5{uj + uiq)] 

Z2 1 

H — ^[S{u - U2q) + S{uj + U2q)] > 

(^^^lmxnn{q,uj). (22) 



1 /T 

nn \u 



This expression is valid for Bose and Fermi superfiuid gases as well as superfiuid "^He, that is, 
all superfiuids described by an order parameter with a phase and amplitude. Equation (122|) 
shows that in the low-frequency two-fluid hydrodynamic region, Wi = Zi/u^ gives the weight 
of the z**^ mode. In contrast, the weight of the i^^ mode in Bragg scattering [see (jl]) and 
f lT7|) ] is given by Bi, where 

B2^^. (23) 

Ul U2 

Comparing (122]) with fll7l) . we see that the hydrodynamic limit (w ^ T) of S{q^uj) involves 
an extra factor of l/u. 

We see that while Imx„„(g,a;) and Imx„„(g, ci;)/^ share the same poles, the relative 
weights of first and second sound are quite different in these two functions, with B2/B1 = 
(W2/Wi){u2/ui) generally being much smaller than W2/W1 due to the smallness of M2/W1 
(see Fig. p. 

The two-fluid expression in fl22p is known from earlier discussions of Brillouin scattering 
in superfiuid ^He (see section HVl) . In contrast to high-frequency inelastic neutron scattering, 
Brillouin scattering is carried out with light in the visible part of the spectrum (w ^ T). As 
pointed out above [see (l22l) ]. this means that Brillouin scattering measures the imaginary 
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part of the density response function divided by the frequency uj. In studies of superfluid ^He, 
this factor of l/cu is crucial to probing second sound since, although only weakly coupled to 
density, the speed of second sound is small and hence, its weight in Brillouin light scattering 
can be significant. 

As noted in section [ITl a similar situation arises in studies of pulse propagation in dilute 
gases, albeit for very different reasons. In this situation [see ([6])], the amplitude of the 
density pulse induced by a localized (in time and space) density perturbation is proportional 
to Imx„„(g, The factor of l/uj here results not from the detailed balance factor, but 

from the Fourier transform of the step-function used to model the time dependence of the 
perturbing laser beams. Arahata and Nikuni [l^ calculated the effect of producing a density 
fluctuation by a sudden perturbation in a superfluid Fermi gas at unitarity. The final result 
was two pulses moving with the speeds of first and second sound, with relative amplitudes 
given by Wi as defined in (!20|) . 



IV. SUPERFLUIDS WITH SMALL THERMAL EXPANSION 

In section UTTI we showed that the relative weights in S{q,uj) and Imx„„(g,ci;) of first and 
second sound are given by Wi = Zijuf and Bi = Zi/ui, respectively, where Zi is defined in 
f llSp . These results are valid for any type of superfluid, including weakly interacting atomic 
Bose gases and also strongly interacting Fermi gases. In this section, we concentrate on 
strongly interacting Fermi gases where, as discussed above, the thermal expansion is rela- 
tively small, meaning that density and temperature fluctuations are not strongly coupled. 
Using the exact relations in (fT4|) . in this section we will show that to a very good approx- 
imation, W2/W1 is given by f l30|) in this situation. In section |Vl we give a brief discussion 
of the dynamic structure factor in a dilute Bose gas, where the thermal expansion can be 
much larger when the s-wave scattering length is small. 



Solving the coupled equations in f|T^ . one obtains the expansions |l4j. Il5| 



ui = vt[l + {^-l)x + ---], M^ = -[l-(7-l)x + ---], (24) 

in terms of the parameter defined by a; = v'^/'jVg. These results show that to lowest order in 
the parameter x(7 — 1) = xeLp, the first and second sound velocities are well approximated 
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o 0.5 h 



0.4 - 
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first sound 




second sound 




TIT 



c 



FIG. 1: Speeds of first and second sound in a unitary superfluid Fermi gas. The dashed 
the approximations ci = Vs and C2 = v/ ^/7, as given by the leading terms in (p^ (from [8 



mes are 



by 



2 2 _ _ T^^o P^o 



(25) 

7 Cp PnO 

In figure [H we show the temperature dependence of ui and U2 and compare these results 
with the leading order expressions Ci and C2 given by (!25l) (shown by the dashed lines). 
The thermodynamic functions needed to obtain these sound speeds are calculated using the 
microscopic theory of Nozieres and Schmitt-Rink (NSR) 25|, l26|]. We see that Ci and C2 are 
extremely good approximations at all temperatures. At temperatures T > 0.4Tc, x = c^/cf 
is very small, even though eLp ~ C^(l)- At lower temperatures, while x is no longer small 
(x = 1/3 at T = 0), eLp becomes extremely small. Thus in both limits, we find that the 
correction term xeLp in flM|) is negligible. 

The fact that -Ui ~ ci and U2 ~ C2 in a fluid with small thermal expansion such as the uni- 
tary Fermi gas means that first and second sound propagate at essentially the same velocities 
as pure density and temperature oscillations, as would arise when eLp = 0. Surprisingly, 
this does not mean that first and second sound are uncoupled density and temperature 
oscillations. We discuss this in detail in [Bl 
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Two-fluid hydrodynamics leads to the coupled equations 



m 



u^5p = 5P, u^6s = sl^5T. (26) 

PnO 

As noted above (see flgure[T]), the speeds of flrst and second sound are well approximated by 
Ci and C2 in ( 125|) . Physically, this means that to leading order, the pressure fluctuations in 
flrst sound are at constant entropy per unit mass s and hence, 6P ^ {dP/dp)sSp. Similarly, 
the temperature fluctuations in second sound are to leading order at constant pressure and 
hence, 6T ~ {dT/ds)p6s. Using these in f l2B]) leads to the leading order expressions for ui 
and U2 given in f l24|) . 

We remark that, qualitatively, our results in flgure[T]for the temperature dependence of 

Ui and U2 are in agreement with Arahata and Nikuni UM, who also based their work on 

I I 

NSR thermodynamics. The differences with [10| at high and low temperatures are easily 
understood. Direct numerical calculations based on NSR are difficult to do accurately for 
T < 0.4Tc because of the extremely small value of the normal fluid density. Our results in 
flgure[T]in this region are based on the assumption that Goldstone phonons are the dominant 
thermal excitations at low T at unitarity, which allows us to obtain the low temperature 
values of Ui and U2 analytically. In the opposite limit of high temperatures close to Tc, our 
calculations show that W2 steadily increases, in contrast with that of p^. For T close to Tc, 



it is crucial to use a superfluid density with the correct critical behaviour 27|, removing the 



spurious flrst-order behaviour predicted by NSR 



(In fact, this behaviour is not unique 



to NSR and is symptomatic of any theory that treats phase fluctuations in a perturbative 



way. See [29| for further discussion.) 

Using the values of Ui,U2, and v (see flgure[T]), we can calculate Z2,W2 and Wi. The 
results are shown in flgure O The relative weight of second and flrst sound in S{q, u) is 
given by 



W2 ^ Z2uf ^ v' - uj 
Using the exact two-fluid expressions for ui and U2 gives the ratio W2/W1 shown in flgure [31 



For comparison, we also plot the Landau-Placzek ratio eLp = 7 — 1, calculated direct 



:tly : 



rom 



the same thermodynamic functions used to obtain the sound speeds Ui, U2, and v 

In flgurelU we plot S{q, u) given by f l22p . This shows the flrst and second sound resonances 
at a series of temperatures. The results shown in flgures |2]|1] show the remarkable feature 
that, even though the maximum value of Z2 is ~ 0.05 at T ~ 0.8Tc, the relative weights of 
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FIG. 2: The temperature dependence of the normahzed amphtudes Zj, Sj, and Wi of first (i = 1) 
and second sound [i = 2) in a Fermi gas at unitarity. 




r/r 



FIG. 3: Comparison between the ratio of the second and first sound amplitudes {W2/W1) and the 
Landau-Placzek ratio eLP = 7 — 1 in a unitary Fermi gas. 
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co/(^Vp) 

FIG. 4: The two-fluid dynamic structure factor given by (j22p at unitarity. The plot shows the 
first and second sound resonances as a function of the frequency oj/qvp (where we take q = O.l/cj? 
and vp is the Fermi velocity) for a series of temperatures from T = to T = Tc, in steps of O.OSTc 
(each temperature is offset). For clarity, the delta functions in (|22p are plotted as Lorentzians with 
a width A = O.Olgui?. 



first and second sound in S{q,u) can be comparable in the temperature region T > O.STc. 
The relative weight of second sound is instead much smaller in Imx„„ (see figure E]) for 
reasons we have discussed. 

Figure El also shows that W2/W1 is quite well approximated by eLp. It is useful to derive 
this result more analytically. Using (fT4|) in (127|) . we obtain after some algebra 

^-^f4-lV- (28) 

n 

This expression was first given by Vinen [30]. Using the expansion for u\ in (JMll . one can 
reduce ^ to 

^ = eLp(l + 7X + ■■■)'■ (29) 
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FIG. 5: The imaginary part of the two-fluid density response function at unitarity. In contrast to 
the results shown in figured where the relative weights of second and first sound are determined 
by W2/W1, in Bragg scattering they are determined by the smaller ratio B2/B1 [see (123p ]. As a 
result, second sound has a much smaller weight in Bragg scattering that first sound. 



At high temperatures, where x = C2/ c\ <^ 1 (see figure [T]), we have 

W2 



eLp[l + 2(l + eLp)a; + 



clp = 7 - 1- 



(30) 



This result was first obtained almost 40 years ago pjj, 1^, 31 1 in the context of Brillouin 
light scattering [l3| in superfiuid ^He near the critical region close to Tc. For T < 0.4Tc, one 
can see from figure [1] that we can no longer assume x ^ 1 and hence expand fl29|) . 

As noted above, the result in (!30|) was obtained in the classic literature on superfiuid 
^He. The Landau-Placzek ratio eLp is extremely small (10~^ — 10"'^) in superfiuid ^He and 
thus second sound has generally a very small amplitude in the dynamic structure factor. 
However, under pressure and close to Tc, the magnitude of eLp in liquid "^He can be signif- 
icantly increased (to about eLp ~ 0.2) such that a second sound doublet in S{q,u) can be 



measured using Brillouin light scattering techniques [13[ . Such studies Q, 
the correctness of W2/W1 = eLp in liquid "^He at a pressure of P = 25 atmospheres in the 



17( 1 have verified 
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temperature region (Tc — T) ~ 10~^ — 10~^K. 

The Landau-Placzek ratio eLp = 7 — ! was first derived in 1934 [32] to describe the relative 
weight of the thermal diffusion central peak vs. sound waves which are exhibited by S{q, u) 
describing the hydrodynamics of a normal fluid. Thus, fl5CT]) generalizes the Landau-Placzek 



result 



32| to the case of two-fluid hydrodynamics. In a normal liquid above T^, the 



dynamic structure factor given by (122!) reduces to 24 1 

S{q, ^) = ^ I ^ [-^(^ - ^'«-^) + ^(^ + ^sg)] + ^^^(^) I • (31) 

Thus, above T^, the second sound doublet in fl22l) collapses into a zero frequency central peak, 
with relative weight 7 — 1 compared to a sound wave at u = Vgq. When one includes damping 
due to transport processes, the central peak in f lM]) broadens, describing the thermal diffusion 
relaxation mode of width DtQ^ (where is the thermal diff usivity) . The behaviour of 
S{q, u) as one goes from below Tc to above Tc is smooth but complicated due to various 
sing ularities in the thermodynamic and transport coefficients (for further discussion, see 

m 



V. DYNAMIC STRUCTURE FACTOR OF DILUTE BOSE GASES 

In this section, we discuss first and second sound in a uniform dilute Bose-condensed gas. 
The two- fluid hydrodynamics is again described by the results given in section IIIII However, 
we need to distinguish between the cases of weakly and strongly interacting Bose gases. The 
latter case can be easily realized by considering a molecular Bose condensate on the BEG 
side of the BCS-BEC crossover in a two-component Fermi gas. In such a molecular Bose 
gas with density ub = n/2 and molecular mass rriB = 2m, the s-wave molecular scattering 

length is given by ~ O.Gap, where ap is the atomic s-wave scattering length between the 

j 

two Fermi components [33|. We will show that the behaviour of first and second sound in 
strongly interacting Bose gases is very similar to that in Fermi gas superfluids near unitarity, 
but very different from that in weakly interacting Bose gases. 

In a dilute, weakly interacting Bose gas, the thermodynamic quantities that enter the two- 
fluid equations can be solved analytically in powers of the interaction strength g = AnaB/mB- 
To first order in g, the speeds of first and second sound are given by (see, for instance. 
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chapter 15 in 



0; 



2 5T 5'5/2(l) 25(no 2 5"^co 



""1 = ^;^ — m + T;^' ""2 = -— • (32) 

Here, 5'„(z) = Xli^i-^V^" is the usual Bose-Einstein function with fugacity z. In the hmit 
of small g, the normal fluid density p„o reduces to the density mfiQ = p — mrico of atoms 
thermally excited out of the condensate. The superfluid density pso is likewise given by the 
density mricQ of condensate atoms. To lowest order in Uco is the condensate density of the 
ideal Bose gas. 

In figure El we show the temperature dependence of ui and U2, satisfying ( fT4l) . in the case 
of weak interactions {nsa^ = 10~^). The thermodynamic functions used to determine these 
sound velocities have been calculated using the Hartree-Fock-Popov (HFP) microscopic 
model for the thermal excitations [3]. The velocities are normalized in terms of a Fermi 
velocity defined by vp = {Gn'^nBY^^ /m (equivalent to the Fermi velocity of a two component 
Fermi gas through the BCS-BEC crossover). We also show for comparison the values of ui 
and U2 given by the lowest order solutions described by (1321) . One sees that these uncoupled 
modes of oscillations are a good first approximation to the full solutions of the two-fluid 
equations in a weakly interacting BEG at all temperatures outside a small temperature 
region at low temperatures where first and second sound hybridize, exhibiting an avoided 
crossing. At this avoided crossing, the natures of first and second sound are interchanged. 
As T — 0, the first sound velocity coincides with the phonon velocity, as is the case with 
strongly interacting Fermi gases and superfluid '^He. For a dilute gas, this phonon velocity 
is the Bogoliubov-Popov phonon velocity ^/grico/mB, and corresponds to an oscillation 
of the condensate component. Above the hybridization point, to a good approximation, 
second sound propagates with this velocity. Consequently, second sound in a dilute Bose 
gas (above the hybridization temperature) is essentially an oscillation of the condensate at 
all temperatures, with a static thermal cloud. In contrast, first sound describes an oscillation 
of the thermal cloud in the presence of a stationary condensate. 

We note that HFP approximation we have used to calculate thermodynamic quantities 
has the well-known problem near Tc of predicting a spurious first-order transition (see, for 



instance, 



34|). As touched on in section llVt this is actually a very general problem, and 



arises in any theory that treats bosonic phase fluctuations of the order parameter as a small 
parameter, including the HFP and NSR theories. Of course, this includes all perturbation 
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BEC 



FIG. 6: Speeds of first and second sound in a weakly interacting molecular Bose gas with gas 
parameter nga'g = 10^^. The blue dashed and red dotted lines are respectively the approximate 
speeds ui and U2 given by ([32|) . The inset shows the LP ratio elp (red dashed line) and the weight 
ratio W2/W1 (black sohd hue). 



ations can 



35| avoids 



theories that go beyond mean-field and only ab initio Quantum Monte-Carlo calcu 
avoid this problem. (A competing fluctuation theory of the BCS-BEC crossover 
this problem by ignoring phase fluctuations.) The unphysical flrst-order transition leads to 
the result that the second sound velocity does not vanish at T^. as it should. In contrast to 
figure [H we have not corrected the spurious behaviour of the superfluid (condensate) density 
close to Tc in the results shown in figures [6] and [71 These figures also show the expected 
feature that the error becomes larger with increasing interaction strength. In spite of the 
shortcomings of HFP theory in the critical region close to T^, we emphasize that it gives a 
good estimate of thermodynamic quantities, outside the critical region. All numerical results 
shown in figures ([6])- ([8]) are obtained using this theory (apart from the uncoupled solutions 
shown in figure [6l which use the ideal gas expression for the condensate density). 

In figure El the inset shows the much larger value of the LP ratio eLp (compared to a 
Fermi superfiuid-see figure and the ratio W2/W1 of the amplitudes of second and first 
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sound in S{q,co) as a function of the temperature. We note that ratio W2/W1 is still given 
by the expressions in (|27|) and fl28l) . However, now the parameter eLp is no longer small 
1) and thus the expansion based on (^^, which leads to (12^ and (|5Ui) . is no longer 
valid. We also explicitly note that W2/W1 is not equal to in either the weak or strong 
coupling limits. 

First and second sound in a dilute, weakly interacting Bose gas involve weakly coupled 
oscillations of the thermal cloud and condensate, respectively. Both modes contribute to the 
density response function. Consequently, both first and second sound are sensitive to density 
probes. This fact has been made use of in a recent experiment involving trapped atomic Bose 
gases by Meppelink et al. [4|]. In this study, a density pulse was generated (as described in 
section [IT]) and the velocity of the second sound pulse was measured above the hybridization 
point by tracking the propagation of the density pulse. The inset in figure [6] shows that 
the amplitude of second sound in S'(g, u) and in density pulse propagation experiments 
is twenty times as strong as first sound for all temperatures above the crossing point at 
T ~ O.lTc. This result may explain the absence of a first sound pulse in the data of jj], 
which corresponded to n^a^ ~ 5 x 10~^. 

In figure [71 we plot first and second sound for a more strongly interacting Bose gas (i.e., 
with a larger value oi g.) This situation is close to the unitarity Fermi gas results shown 
in figure [1] 36|] . In particular, we see that the velocities of first and second sound never 
become degenerate (cross) and thus the hybridization shown in figure [6] does not occur. In 
this strongly interacting Bose superfluid, we also show the approximate velocities of ui and 
U2 given by (1251) . These approximate values give a reasonable first order approximation to 
the full solutions of the first and second sound velocities, rather than the expression in ([32]) 
that applies to a weakly interacting Bose gas as shown in figure [61 The inset shows the ratio 
of W2/W1 vs. T, as computed from (127]) . In this case, the amplitude of second sound in 
S (g, u) is only 2-3 times larger than first sound for T > 0.5Tc. 

In figure [HI we compare the values of Wi/{Wi + W2) and W2/{Wi + W2) for a weakly 
interacting Bose gas {nBa\ = 10^^) and a strongly interacting Bose gas {ubcl^ = 0.1). One 
sees that the latter case is qualitatively similar to the result in a Fermi superfluid gas at 
unitarity given in figure [2l In contrast, in a weakly interacting Bose gas, for the temperature 
range above the hybridization point, the amplitude of second sound is much larger than first 
sound in both S{q, co) and density pulse propagation experiments. 
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TIT 

BEC 

FIG. 7: Speeds of first and second sound in a strongly interacting molecular Bose gas with gas 
parameter n^a'g = 10^^. The blue dashed and red dotted lines are respectively the approximate 
speeds ci = fg and C2 = v/y/^ given by ((25]) . The inset shows the LP ratio eLP (dashed line) 
and the weight ratio W2/W1 (black line), to be compared with the result in figure [3] for a strongly 
interacting Fermi gas at unitarity. 



VI. BRAGG SCATTERING IN A TRAPPED ATOMIC GAS 

For illustration, we now consider the response of a trapped Fermi gas at unitarity to two- 
photon Bragg scattering due to first and second sound. We consider an experimental setup 
similar to that used in 37| , where a pair of laser beams are applied to the gas such that the 
region where they overlap is confined to the centre of the gas. As a first approximation to 
this configuration, we assume that the overlap region between the two beams is isotropic, 
with a Gaussian profile: 

6V{v, t) = -^e-^'/2a2 ^^g^q . r - u;t) . (33) 
v27rcr 



Here, a is the width of the Bragg beams and q and u are the wavevector and frequency 
difference between the two beams. 

For a Bragg pulse of short duration (r ^ where ojq is the trap frequency), the 
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FIG. 8: Temperature dependence of the normalized amplitudes Wi/{Wi + W2) (solid lines) and 
W2/{Wi + W2) (dashed lines) for a weakly interacting Bose gas (n^a^ = 10^^) and a strongly 
interacting Bose gas (n^a^ = 0.1). 

momentum transferred to the gas is related to the imaginary part of the density response 
function [see (jlj)], through a convolution of the atomic density: 



In this expression, Imx„„(g, r) is the imaginary part of the density response function, 
evaluated within a local density approximation (LDA). The imaginary part of the two- fluid 
density response function is given by ffT7|) and should be evaluated at the local value of the 
density. The LDA is justified for wave vectors much larger than the inverse of the size of 
the trapped gas, so that the system can be considered, locally, as a uniform body. 

Results for the momentum transferred given by flM|) . as a function of u/qvp, are shown 
in figure [3 for T = 0.75Tc and choosing a = 0.5Rtf, where Rtf = \/h/muj{24Ny^^ is 
the Thomas-Fermi radius of a noninteracting two-component Fermi gas. Although small, 
the second sound peak is well separated from the first sound peak, revealing that two- 
photon Bragg spectroscopy might become a practical way to measure second sound. Above 
and below this temperature, the contribution of second sound to the momentum transfer 




(34) 
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FIG. 9: Momentum transferred to a trapped Fermi gas at unitarity by two photon Bragg scattering 
due to first and second sound at T ~ 0.75Tc. Solid (blue) line shows the momentum transferred due 
to first sound as a function of the beam detuning a; divided by the relative wave vector q= | qi — q2 1 ■ 
vp \& the Fermi velocity. Dashed (red) line shows the momentum transferred due to second sound. 



becomes significantly smaller, however. 



VII. CONCLUSIONS 



In this paper, we have given a detailed analysis of the density response function \rin 
and the related dynamic structure factor S{q, co) for uniform superfluids in the two-fluid 
hydrodynamic limit. We have put special emphasis on the second sound contribution to 
S{q,u) in the case of a strongly interacting Fermi gas superfiuid (at unitarity). While 
many features of the dynamic structure factor we discuss are known in the superfiuid "^He 
literature [is], Q , several new aspects arise in the case of ultracold superfiuid gases. 

First and second sound both appear as resonances in the density response function but 
their relative weights vary in different experimental probes. Two-photon Bragg scattering, a 
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well-developed tool in ultracold gases used to study density fluctuations (so far, mainly in the 



collisionless region 



20 



38 



39|), is directly proportional to Imx„„(g, u). In contrast, the low- 



frequency limit of the dynamic structure factor S{q,u) = — Imx„„(g, u;)/7rn[l — exp{—f3u)] 
measured in Brillouin light scattering, and also the amplitude of density pulses that can 
be generated in dilute gases [lO|, are proportional to Imx„„(g, w)/^. This extra factor of 
l/u leads to a large enhancement of the weight of low frequency second sound in S{q,u) 
compared to flrst sound. 

The relative spectral weights of flrst and second sound in x„„(g, w) are very dependent 
on how strongly temperature fluctuations couple to density fluctuations. For S{q,u), this 
coupling is conveniently parametrized by the value of the dimensionless Landau-Placzek 
ratio eLp = 7 — 1, where 7, deflned in f lTB]) . is the ratio of the speciflc heats (per units mass) 
at constant pressure and volume. At flnite temperatures, the relative weight of flrst and 
second sound (W2/W1) is equal to eLp in both unitary Fermi gases and superfluid ^He. The 
key difference between superfluid ''He and a unitary Fermi superfluid lies in the value of eLp. 
Apart from a small region very close to the transition temperature Tc, one has eLp ^ 1 in 
superfluid ''He Q, In contrast, for T > 0.7Tc, one flnds that eLp is of order unity in 
Fermi superfluids (see flgures[3]and|l]) and hence the second sound resonance has appreciable 
weight in S{q, u). 



Density pulse experiments as suggested in 



10| seem to be a promising way of detecting 



both flrst and second sound pulses in Fermi gases, although more work is needed to clarify 
the nature of flrst and second sound in cylindrical geometries where such pulses would be 
generated. (See, for instance, recent work on this problem in 40|].) The experimental results 
reported in [2^] , where a blue-detuned laser was used to generate a density pulse in a Fermi 
gas close to unitarity, revealed no clear signs of a second sound pulse. This experiment was 
carried out at very low temperatures, however, where second sound will not have appreciable 
weight in S{q,co) (see flgures [3] and H]) . Alternately, even though the weight of second sound 
is comparatively small in Bragg scattering, Bragg scattering has the advantage that it can 
be used to probe in a very precise way the properties of uniform superfluid Fermi gases js^] , 
including the velocity of second sound. 

In section |Vl we compare the two-fluid hydrodynamics in a uniform, dilute, weakly in- 
teracting Bose condensed gas (flgure [6]) with a strongly interacting Fermi superfluid gas. 
In such Bose superfluids, the thermal expansion (and hence, eLp) is much larger, with the 
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result that second sound is the dominant excitation in S{q, u) and mainly involves a pure 
oscillation of the condensate in the presence of a static thermal component js, 4|. However, 
first and second sound in a strongly interacting Bose gas are similar to those in a Fermi gas 
near unitarity (see figure [7]). 
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Appendix A: Phonon thermodynamics at low temperatures 



At low temperatures, Goldstone phonons determine the thermodynamics of both super- 
fluid "^He and Fermi gases. For phonons with velocity c, the free energy F and normal fluid 
density pno are given by [4l| (Recall that we have set h = = 1) 



and 



90c3 



45c5 



Using flAip . it is straightforward to show that 

1 fdF 

So 



mN \dT 



45pc3 ' 



(Al) 

Pno = (A2) 

(A3) 

-Po+'-^ 1 + — (^1 , (A4) 

c. = — ^. (A5) 

In arriving at these expressions, the temperature dependence of the Goldstone phonon ve- 
locity c has been ignored. Pq is the pressure in the ground state. These results can be 
combined to give 



dv 



T,N 



and 
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and 



V = s 
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(A7) 
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Combining the above results, the Landau-Placzek ratio at low temperatures is given by 

2 



CLP 



1 + 



3p f dc\ 



and hence, recalling that x = v'^/'-fVg and 7 = Vg/v^, 
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3p /^5c\ 
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This result is valid for both the superfluid unitary Fermi gas and superfluid ^He. Further 
simplifications are not possible without ab-initio or experimental information about c. At 
low temperatures, however, where p^o — Po and Vg — c (see figured]), (lA9p reduces to 



PsO 

( 

PnO 



3p /dc\ 

C \dpjr 



(AlO) 



Equations (IA6p and (]A7p also give us the result that 7 — )■ 1 and x = v^l^v\ — )■ 1/3 as 
T — 7- 0. Finally, for a Fermi gas at unitarity, {p/ c){dc/ dp)T = 1/3 and flAlOp is further 
reduced to the result given in fIBlOp . In superfiuid ^He, the Griineisen constant is almost a 
factor of 9 larger, {p/c){dc/dp)T ^ 2.84 |42|. 



Appendix B: Superfluid and normal fluid velocity flelds in second sound 



As Landau [2| originally noted, when the thermal expansion coefficient can be neglected 
(i.e., eLp = 0), one can show that second sound is a pure temperature oscillation, with no 
associated density fiuctuations. This last feature follows from the fact that when eLp = 0, 
second sound involves no mass current (j = Pno^n + Pso^s = 0) and thus the continuity 
equation requires that 6p = 0. Likewise, first sound is a pure density oscillation, and involves 
no fiuctuations in the entropy per unit mass s. The non-dissipative two-fiuid equations give 
the following equation for this quantity: 

^ = — — V-(v, -v„). Bl 
dt Po 

This shows that for pure adiabatic motion {6s = 6T = 0) such as first sound, the superfiuid 
and normal fiuid velocities are identical: = v„. Putting these results together, the velocity 
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fields for first and second sound when eLp = are given by 



PsOV. 



(2) 



-' (2) 
PnOVn 



(B2) 



where, for instance, denotes the superfluid velocity field associated with first sound 
(along the direction of propagation q). This Appendix discusses how these velocity ratios 
change when eLp is non-zero. Surprisingly, we find that changes are quite substantial at all 
temperatures, including the low-T region where eLp ^ 1. 

As we discuss in the present paper and in sl, the Landau-Placzek ratio eLp increases with 
temperature and in Fermi gases near unitarity, eLp is of order unity at high temperatures 
(see figure [3]). As we have shown, this is large enough that second sound has significant 
weight in the dynamic structure factor. 

To understand the effect that a non-zero value of eLp has on the superfluid and normal 
fluid v„ velocity fields, one can solve the two-fluid equations for the values of these fields 
associated with first and second sound u = uq, where u = ui oi U2- Defining the variables 
j = PsqVs + Pno^n and w = pso(vs — v„), the plane- wave solutions of the linearized two-fluid 
equations are 



Here, 
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has units of velocity squared. Recall that x is defined below fl24|) . Standard thermodynamic 
identities have been used to introduce eLp = 7 — 1. As the structure of (IBSp makes clear, 6^ 
determines the coupling between density (oc j) and entropy per unit mass (oc w) fluctuations. 
[The latter is given by (IBip .] 

Using ( IB3p . it is straightforward to derive an expression for the ratio of the velocity fields 
associated with first and second sound: 
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Substituting the expressions for the speeds of first and second sound given by IHM into ( IBSp . 
one finds 

(1) 2 1,2 

(B6) 
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and 

P„ot'f " .1(1-.) -6= ■ * ' 

In section HVt we showed that ci and C2 in fl25|) . which are expansions in the small parameter 
eLpx, are excellent approximations to the speeds of first and second sound at all temperatures 
(see figure [I]). We thus expect ( 1B6I) and ( 1B7I) to be similarly good approximations for the 
velocity fields at all temperatures. Using flB4p to remove the dependence of these expressions 
on Vg, they reduce to 

(1) 1 - Jf^e^px 

- V (B8) 

and 



PnOVlr (1 - a;) - ^CLPX 

When eLp is set to zero, we recover the results in (]B2p . describing pure in-phase and out- 
of-phase oscillations of the superfiuid and normal fiuid components. However, we note that 
even though eLp vanishes as T — )■ 0, so does the normal fiuid density pno, and it is not 



obvious that a/ Pso^Lpa^/ Pno is small at low temperatures. Conversely, at high temperatures. 



T — Tc, the factor of a/p^io7Pso multiplying y/ei^-px in (IBSP and flBOp will enhance the small 
parameter ^eLpx. It is thus crucial to know the temperature dependencies of Pso/Pno and 
eLpx in both limits. 

At low temperatures, the thermodynamics is dominated by Goldstone phonons. Using 
standard results in the phonon regime, one can show (see [X]) that in the T — limit. 



^eLPX ^ I (BIO) 

PnO J 

for a unitary Fermi gas. This non-zero limiting value as T — )■ arises from the fact that both 
eLpx and p„o vanish as T^. For a Fermi gas at unitarity in the T — > limit, using flBlOp and 
the fact that x = 1/3 (see dHS]) and give 

(1) Q (2) 

— -7' — ^--2, (Bll) 

even though eLp vanishes as T — )■ (as shown in figure [3]). If we had simply set eLp = 0, as 
often done in the superfiuid hterature, we would instead obtain the results given by ( ]B2p . 

In figure [10], we plot the velocity ratios in (1B8I) and (1B9|) for a superfiuid unitary Fermi 
gas. It is clear that the superfiuid and normal fiuid velocities which are involved in first 
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FIG. 10: Superfluid and normal fluid velocity fields associated with first and second sound in a 
unitary Fermi superfluid. The solid (black) line shows the ratio v^g'' /vn'^ of the superfluid and 
normal fluid velocity fields involved with first sound, ()B8|) . The dashed (red) line shows the ratio 
PsQVs / Pnovii of the superfluid and normal fluid currents associated with second sound, ()B9p . If 
we set eLP = 0, these ratios would be 1 and —1, respectively, describing the velocity fields for pure 
density and temperature oscillations. 



and second sound are significantly different from those given in flB2l) . This is a reflection of 
the fact that flrst and second sound are not pure uncoupled density and temperature waves, 
respectively, which would be described by ( IB2I) and correspond to eLp = (and p„o 7^ 0). 
As we have shown in Q], such in-phase (v^ = v„) and out-of-phase zero current (j = 0) 
solutions correspond to flrst and second sound velocities given by 



Ml 



ul 



T^l2l=v^. (B12) 

|-| PnO 

As we show in flgure 2 of [8|, these sound velocities are in good agreement with a full 
calculation including the coupling associated with eLp. In this Appendix and section HVl we 
have carried out a careful analysis based on ex pan ding to flrst order in the parameter ei^px, 



which is small (^ 1) at all temperatures IJ, ll5|. The results in flgure [TO] show that the 
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main effect of working with 7 not equal to unity is simply that the second sound speed U2 is 
given by vj ^ as in fl25|) . instead of v in flB12|) . In contrast, the first sound speed is always 
very well approximated by fg, even when 7 deviates from unity. 

Initially, it seems surprising that the superfiuid and normal fluid velocity fields shown in 
figure [in] are so different from those in (1B2I) . since the first and second sound velocities are 
fairly well approximated by ( ]B12p . (See figure 2 of sj].) This arises because of two features 
evident in ( IBSP and flB9|) . First of all, we note that the first and second sound speeds in 
( 12^ involve corrections of order eLpX <^ 1, while the related corrections are much larger in 
( 1B9P and ( IBllI) since they enter as ^e\spx. The second, and more important, reason is that 
the corrections involve the factors ■>/ Pso / Pno and a/ Pno/Pso, which are very large at T = 
and T = Tc, respectively. 

Substantial deviations from the ratios given in flB2|) also occur in superfiuid "^He, although 
the magnitude will be different from the case of a Fermi superfiuid at unitarity shown in 



fig 
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ure 



3l| 



TOl This feature was not commented on in the older superfiuid He literature 



which concentrated entirely on the frequency and damping of the first and second 
sound resonances appearing in the dynamic structure factor. The associated velocity fields 
calculated here were not discussed. These velocity fields may be direct experimental interest 
in future studies of two-fiuid hydrodynamics in Fermi superfiuids. 

Calculation shows that the various terms in flB8l) and (IB9|) are not small compared to 
unity. However, if one formally expands these expressions to leading order in (= 



\/pnoeLpa;/p5o), one finds 



(1) (2) 

^ = 1-5, ^ = -1-5, (B13) 

Vn PnQVn 



where the leading correction to (IB2p is found to be 



d = — — = — J— eLpx. (B14) 



PnO PnO V PsO 

This value of S agrees with the expression originally worked out by Lifshitz (see page 71 in 
43|) in the limit where u"^ <^ uf, valid close to Tc. In the case of superfiuid ^He, the fact 
that eLp ^ 1 even close to (it diverges very weakly at T^) means that the first correction 
in ( ]B13|) is indeed small, 5^1. In contrast, for a Fermi gas at unitarity, eLp is of order 
unity or larger near Tc (see figure [3]) and as a result we find S = 0.3 just below Tc [m]- In 
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this case, the expansion in (1B13|1 is not vahd, as can be seen by comparing (IBlSp with the 
results in figure [TOl 
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